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Abstract 

Ordering processes in fcc-alloys with composition A%B (like CU3AU, CU3PCI, CoPt3 
etc.) are investigated by Monte Carlo simulation within a class of lattice models 
based on nearest-neighbor (NN) and second-neighbor (NNN) interactions. Using an 
atom-vacancy exchange algorithm, we study the growth of ordered domains follow- 
ing a temperature quench below the ordering spinodal. For zero NNN-interactions 
we observe an anomalously slow growth of the domain size L(t) ~ t a , where a ~ 1/4 
within our accessible timescales. With increasing NNN-interactions domain growth 
becomes faster and a gradually approaches the value 1/2 as predicted by the con- 
ventional Lifshitz-Allen-Cahn theory. 
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1 Introduction 



Ordering transitions in real metallic alloys normally show kinetic properties 
which are substantially more complex than the standard Lifschitz-Alien-Cahn 
scenario in the case of a scalar, non-conserved order-parameter [1,2]: i) The 
symmetry of the ordered phase can imply a multicomponent order-parameter 
and the appearance of several antiphase domains, ii) There may be differ- 
ent types of domain boundaries providing different forces for curvature-driven 
coarsening, iii) Coupling of non-conserved order-parameter components to the 
(conserved) alloy composition generally leads to compositional changes within 
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domain boundaries. Domain growth hence is accompanied and slowed down 
by interdiffusion processes, iv) Depending on properties of the chemical inter- 
actions, vacancies may enrich in the domain boundaries and thus enhance the 
dynamics when atom-vacancy exchange is the prevailing migration mechanism 
[3,4]. v) Quenching into a two-phase region in the temperature-concentration 
plane releases competing processes of ordering and spinodal decomposition, 
eventually accompanied by the appearance of transient phases, vi) Finally we 
mention effects caused by surfaces or imperfections, which break the lattice 
translational symmetry. 

Despite a large body of literature [5], many questions in this area have re- 
mained open. Our aim here is to study some new qualitative aspects primarily 
in connection with the first two points listed above. For that purpose we con- 
sider fee A 3 B-type alloys (such as Cu 3 Au, Cu 3 Pd etc.) that display an ordered 
LI2- structure. A class of models is considered which includes chemical inter- 
actions on the fcc-lattice up to second neighbors. Our main result will be to 
demonstrate a gradual changeover from an anomalously slow growth with an 
effective exponent a ~ 1/4 to conventional growth, characterized by a = 1/2, 
when the second neighbor interactions are increased from zero. Results will 
be interpreted in terms of properties of low-energy or type-I domain walls, 
whose energy vanishes as the second-neighbor interactions are switched off. 



2 Fee— alloy model with vacancy— driven kinetics 



Consider an fcc-lattice, where each site i can be occupied by an A-atom, a 
S-atom or a vacancy (V). Occupation numbers thus satisfy cf + cf + = 1. 
Their averages over all sites are chosen in accord with stoichiometric A 3 B 
alloys, c A = 3c B . The average vacancy concentration c v <C 1 is taken small 
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Fig. 1. Comparison of domain patterns reflected by ipi at T = 0.5 Tq for a) R = 0; 
t = 8 • 10 3 MCS and c) R = 0.3; t = 3 • 10 2 MCS, while ^3-patterns are shown in 
b) and d) for the same parameters. From the figures type-I and type-II walls can 
clearly be distinguished. 
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enough so that static properties are unaffected by the vacancies. For NN-sites 
% and j interactions are taken as in our previous work [6], V BB = V BB > 0; 
V^ A = Vjj B = —V BB . On the other hand, for NNN-sites i and j we assume 
V? B = V AA = RV BB ; V AB = 0, with < R < 0.5. 

The above model displays a 4-fold degenerate ground state corresponding to 
the ideal Tl 2 -structure. This structure is described by a 4-component order 
parameter [7], \1/ = (ip , ipi, ip 2 , ^3), where ip refers to the ^/^-concentration, 
while tpi, tp2 and 1^3 describe a succession of 5-rich and 5-depleted atomic 
layers along the x-, y- and z-direction, respectively. Ordered domains are 
of the type (^, oc (-1, 1, 1); (1, -1, 1); (1, 1, -1) and (-1, -1 - 1). A 

sign-change of ipi along the x-direction implies a high-energy type-II domain 
wall. Along the y- or ^-direction a sign change of ipi can be effected without 
braking nearest-neighbor bonds, and creates a low-energy type-I wall. Its 
energy solely results from NNN-interactions and therefore is zero if R = 0. 

In our Monte Carlo simulations we employ the atom-vacancy exchange algo- 
rithm. The ordering temperature To and the temperature for spinodal order- 
ing T sp are deduced from simulations as in [8]. At R = we have ksT = 
1.83Vbb/2, while T sp ~ 0.967 T . T and T sp increase nearly linearly with R, 
such that T o (R)/T o (0) ~ 2.38 and T sp (R) /T sp (0) ~ 2.43 at R = 0.5. 

Typical domain patterns observed after a sudden quench from infinite tem- 
perature to a final temperature T < T sp are shown in Figs. la,b for R = and 
in Figs. lc,d for R = 0.3. Evolution times after the quench were chosen such 
that the typical domain sizes (see section 3) nearly agree in both cases. Apart 
from the much faster evolution in the case R = 0.3 the most important obser- 
vation is that type-I walls in Figs. la,b are flat and nearly perfect, whereas in 
Figs. lc,d they show curvature and larger fluctuations. A connection of these 
different behaviors to kinetic properties will be discussed in section 3. 

Using experimental tracer diffusion constants for calibration, our algorithm 
allows us to relate the Monte Carlo time to the physical timescale in kinetic 
processes. Details of this analysis, applied to Cu 3 Au, are described in [6]. As an 
order-of-magnitude approximation we find that 1 MCS roughly corresponds, 
for example, to 0.1s at T = 0.9 To and 10 7 s at T = 0.5 To. We come back to 
this point in section 3 when discussing different temporal regimes in domain 
coarsening processes. 



3 Domain growth results and discussion 

For a quantitative analysis of the structure and growth of domains we intro- 
duce equal-time structure factors S a (k,t) = {\^f a (k, t)\ 2 }. In view of the sign 
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changes in the local order parameters ip a (r,t) across a wall, see section 2, it 
is clear that the linear combinations 



S\\(k, t) = ± (S^k, 0, 0, t) + S 2 (0, k, 0, t) + S 3 (0, 0, k, t)) (1) 

and 

S±(k, t) = — ]T ( S i(°> 9i» 92, t) + S 2 (q u 0, q 2 , t) + S 3 ( qi , q 2 , 0, f)) (2) 

reflect the arrangement of type-II and type-I walls, respectively. In Eq. (2), Nj, 
denotes the number of pairs (q± , q 2 ) that fulfill q\ + q\ = k 2 . Typical distances 
between walls of either type are determined by the first moments ku(t) and 
k±(t) of the structure factors (1) and (2). In addition we study the excess 
energy stored in the domain walls, AE(t) = E(t) — E(oo), where E(oo) is the 
energy after complete equilibration. 

In the following we focus on T < T sp . In Fig. 2a-d Ai?(t)-data are presented 
for different R. At the highest temperature T = 0.9 T , which in all cases is 
slightly below T sp , a shoulder as a remnant of an incubation period in the 
metastable regime is still visible, while a power-law-type decay of AE(t) pre- 
vails for deeper quenches. A similar behavior is found for ku(t), see Fig. 3 and 
also for k±(t) (not shown), which is always larger than k\\{t). The important 
observation is that within the timescales of our simulations these power-laws 




Fig. 2. Evolution of the excess energy AE(t) per lattice site for different tempera- 
tures T and different NNN-inter actions: a) R = 0; b) R = 0.1; c) R = 0.3 and d) 
R = 0.5. Dotted lines have slopes corresponding to a = 1/4 and a = 1/2. 
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change with R. As mentioned already, at R = where type-I walls have zero 
energy, an anomalously slow growth is found, that can be represented by an 
effective exponent a ~ 1/4 within an extended time window up to our largest 
computing times. Upon introducing small second-neighbor interactions with 
R = 0.1 we observe a significantly faster growth. This becomes particularly 
evident by comparing the data for the lowest temperature T = 0.5 To in Figs. 
2a and 3a with those in Figs. 2b and 3b. The slope of these data corresponds to 
exponents still significantly smaller than 1/2. Further increase of R to R = 0.3 
(Figs. 2c and 3c) yields exponents close to the "normal" behavior a = 1/2, per- 
taining here to T < 0.7 T and to a time regime starting even below 10 2 MCS. 
No substantial change occurs when going to R = 0.5, as shown in Figs. 2d 
and 3d. The steeper decay near t ~ 10 4 MCS can be traced back to finite 
size effects, as A;^" 1 becomes of the order of the size L = 128 of the simula- 
tion cell. These findings support the conclusion that the extraordinary slow 
growth in the case of zero NNN-interactions (R = 0) arises from the existence 
of zero-energy, curvatureless domain walls of type I, which are extremely sta- 
ble. Within the time window considered, such behavior appears representative 
for a class of systems with R small, where type-I walls have sufficiently small, 
but non-zero energy. 

Clearly, our results especially for R < 0.1 do not allow us to draw any con- 
clusion as to the exact asymptotics for quantities displayed in Figs. 2 and 3. 




Fig. 3. Evolution of the first moment (t) of the structure factor 5|| (k, t) for different 
temperatures T and different NNN-interactions: a) R = 0; b) R = 0.1; c) R = 0.3 
and d) R = 0.5. 
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However, from the discussion at the end of section 2 it appears that simula- 
tions up to 10 4 MCS already exhaust the timescales relevant to experiments 
for quenches sufficiently below the spinodal. 

From the discussion of Fig. 1 it is clear that structure factors S a for a given 
a are affected by a 1-dimensional (1-d) array of type-II walls parallel to the 
cv-direction and a 2-d array of type-I walls perpendicular to it. Hence we 
expect the structure factors S\\ and S± to obey scaling laws for 1-d and 2-d 
systems, respectively [9]. This is verified in Fig. 4 showing master curves for 
Su and S± when scaled according to Porod's law in one and two dimensions. 
This plot, applying to R — 0.5, is very similar to analogous plots in [6] for 
R = 0. 
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Fig. 4. Structure factors S\\(k,i) a) and S±(k,t) b) scaled according to Porod's law 
in d = 1 and d = 2, respectively, at T = 0.5 Tq for R = 0.5. 
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